Lensing effects in inhomogeneous cosmological models 
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Concepts developed in the gravitational lensing techniques such as shear, convergence, tangential 
and radial arcs maybe used to see how tenable inhomogeneous models proposed to explain the ac- 
celeration of the universe models are. We study the widely discussed LTB cosmological models. It 
turns out that for the observer sitting at origin of a global LTB solution the shear vanishes as in the 
FRW models, while the value of convergence is different which may lead to observable cosmological 
effects. We also consider Swiss-cheese models proposed recently based on LTB with an observer 
sitting in the FRW part. It turns out that they have different behavior as far as the formation of 
radial and tangential arcs are concerned. 



I. INTRODUCTION 

There is an on-going debate in the community wether 
inhomogeneities on scales up to several hundred mega 
parsecs could account, at least partially, for the observed 
dimming of the SNe la [Celerier [1], Sarkar [2]]. Irre- 
spective of successes of these models to explain the dark 
energy, it is desirable to have some independent cosmo- 
logical tests of these models. The obvious effect of the 
inhomogeneities is the light propagation which may differ 
significantly in a clumpy universe compared with a ho- 
mogeneous FRW model. Recently Vanderveld et al [3], 
being interested in the possibility of explaining the su- 
pernova data, have studied light propagation in one of 
the Swiss cheese models recently proposed [4] using weak 
field gravitational lensing techniques. We are interested 
in the potential gravitational lensing effects in inhomo- 
geneous models to sec if there arc crucial effects which 
could constrain or rule out some of the models. 
In this paper we focus on LTB-bascd inhomogeneous 
models. The LTB solution of Einstein equations de- 
scribes an inhomogeneous spherically symmetric dust 
filled cosmological model with a distinguished center of 
symmetry. Having a center of symmetry, thus not re- 
flecting the homogeneity of the universe at large, these 
models are mainly used as toy models to study the lo- 
cal inhomogeneities of the universe. After a short review 
of LTB inhomogeneous solutions, and Swiss-cheese mod- 
els based on it, we look at the general question of how 
different does a LTB model behave relative to a FRW 
model as far as gravitational lensing concepts are con- 
cerned (section III). For an observer at the origin of a 
flat LTB metric, the shear vanishes as in the case of a 
FRW metric. However, the convergence is different from 
that of the FRW due to the r-dependence of metric coeffi- 
cients and the corresponding scale factor. The vanishing 



of the shear is only true for an observer at the origin of 
LTB. The more general and interesting case of an off- 
center observer has to be treated separately. This case, 
however, maybe treated in Swiss cheese models based on 
LTB which is the subject of the section IV, where arcs 
in strong lensing regime are investigated for two differ- 
ent Swiss-cheese models, assuming the observer and the 
source are sitting in the cheese (FRW metric). While 
the onion model [5] turns out to produce neither tangen- 
tial nor radial arc, the model proposed by Marra, Kolb, 
Matarrese and Riotto [4] do produce both radial and tan- 
gential arcs. The models we have considered could be 
looked at as providing different density profiles of a lens, 
to be compared with the profiles published so far and 
mainly based on simulations. The Swiss cheese models, 
e.g., provide us with density profiles which maybe com- 
pared for special parameters to the NFW one which is a 
phenomenological density profile well supported by the 
simulations [6]. 



II. A REVIEW ON INHOMOGENEOUS LTB 
MODELS 

The LTB metric (c = 1) in co- moving and proper time 
coordinates with zero cosmological constant is given by 

ds^ = -dt^ + S^{r, t)dr^ + R^{r, t){d9^ + sin^ 0d(j)^). 

(1) 

Considering dust stress-energy tensor, Einstein's equa- 
tions imply the following constraints: 
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where dot and prime denote partial derivatives with re- 
spect to t and r respectively. Function p{r, t) is energy 
density of the matter. Functions E{r) and M(r) are left 
arbitrary. M{r) can be interpreted as the mass inside 
of co-moving sphere with coordinate radius r. Assuming 
M (r) > 0, M{r) can be chosen as 
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Now, taking this as the definition of the coordinate ra- 
dius, we may write solutions of Einstein's equations de- 
pending on E{r) in the following ways: 
1. For E > the solution is 
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2. For E — wc have 

i?(r,t) = [^GAf(r)]i/3[i_i„(^)]2/3. 

3. Finally for E < the solution is 
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The so-called bang time function t„(r) is an integration 
constant, indicating different singularities defined by i = 
in [7]. 

As we are going to discuss phenomenon far from these 
singularities, we may assume t 3> tn- Therefore, we will 
assume from now on t„ =0. 
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Keeping next to leading terms in u, we obtain from (11) 
and (12) 
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As mentioned before, the relation between parameter u 
and spatial curvature k{r) is given by 



u = —T\/ k{r) 
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Therefore, the small u approximation is valid when 



= 7Tv/%) < 1, 



(15) 



where 



i?2 := = 0.08, 



Substituting u yields, 

R{r, t) = ^7'tV[1 + i?27V2fc(r)]. (16) 



A. Small-u approximation 

Relation between the coordinate r and the parameter 
u in the non-flat LTB cases is not trivial. It has been 
shown that for the special case of E{r) being a trigono- 
metric function of r, small u approximation is even valid 
for enough large r [5]. This, however, is not in general 
the case as may be seen for a polynomial function E{r). 
However, to simplify the calculation, we are going to as- 
sume u to be a small parameter. This approximation, 
which can describe the dynamics even when 5p/p ^ 1, 
allows to solve the Einstein's equations [5,8]. 

Physically is related to spatial curvature E{r)/r'^ 
[5,8]. In addition we assume > 0, leading to the fol- 
lowing equations: 



B. Swiss-cheese model 

The inhomogeneous metrics may be used to model uni- 
verse in different ways. The direct way is to take an 
inhomogeneous metric, say a LTB solution of Einstein 
equation, as the model universe and see the effect of 
lensing in it. One may, however, devise a so-called Swiss- 
cheese model in which the bulk (cheese) is represented by 
a matter-dominated flat homogeneous FRW model and 
the spherically symmetric holes are constructed using a 
specific LTB solution. The holes which represent the 
inhomogeneities are distributed randomly in the bulk, 
so the model is isotropic and homogenous on average. 
The matching of the inhomogeneous holes to the FRW 
bulk must be handled with care [9] . Depending on differ- 
ent types of LTB solutions, one may construct different 



Swiss-cheese models. 

Biswas et. al [8] study a Swiss-cheese model in which 
the holes arc represented by a LTB metric in the small 
u approximation regime (section II A) . They choose Mq 
in such a way that the coordinate density, Mq, coincides 
with the average density (po) at present time to- 
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The matching conditions imply [10] 
k'{L) = 0, 

where L is the comoving radius of the hole and prime 
means derivative with respect to r. Using the above nor- 
malization (17) for Qf^ <^ 1 (f2fc is curvature abundance 
of the homogeneous universe at the same time)we arrive 
at 
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k{L) = —n^. 

In order to be consistent with CMB we choose 
k{L) = 0. 

According to [10], continuity at the origin implies another 
constraint on curvature: 

k'{0) = 

Recently Marra et.al [4] defined a different Swiss-cheese 
model in which arbitrary number of spherical holes with 
different size and density profile are distributed in the 
cheese. The cheese evolves as FRW while the holes evolve 
differently. At the boundary of the holes, as a conse- 
quence of the boundary conditions, the average mass den- 
sity, defined by p = j^f^^f.y.i Jq p{r,t)R'^R'dr, coincides to 
the FRW density, and E(r) has to go to zero. As far as 
local physics is concerned, the hole has no effect on the 
observer outside it. 



III. LTB UNIVERSE - OBSERVER AT THE 
ORIGIN 

Let us first study the lensing effect due to the global 
inhomogeneity in a flat LTB solut ion (i?(r) = 0): 



ds^ = dt^ - R'{t,rfdr^ - R{t,rfdn^ 



(18) 



We do not assume any single or multiple lens but are in- 
terested in the global effect of bending of light rays due 



to the LTB inhomogeneities. For simplicity, we place the 
observer at the event O (center of the inhomogeneous re- 
gion) with 4- velocity , u'^Uao = 1. Choosing the affine 
parameter of the rays. A, at O such that (1) A = at 
origin, (2) A increases to the past and (3) fcau" \o~ ~1, 
then fc" = ^jy^ is past directed. Using the dimension- 
less fc", the corresponding wave vector is then defined by 
— ^fc", where uJq is the frequency of the wave measured 
by the observer at O. 

Let 7o be a ray and on 70 be the result of the par- 
allel propagated four velocity of the observer, u". The 
orthonormal bases along 70 on the lens plane are and 
E2 . The deviation vectors of the beam centered on 70 
can then be written as = —S,iEf—£_2E2—£,[)k", where 
^1 and ^2 are called the screen components of the corre- 
sponding separation vector of two neighboring light rays 
[11]. For the above metric these vectors are derived as 
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where z is the red-shift of the source (O) defined as 
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The evolution of ^ = (^1,^2) is given by the following 
equation of geodesic deviation [12]: 



|(A)=T(A)C(A), 



(22) 



where T is the optical tidal matrix describing the influ- 
ence of space-time curvature on the propagation of light: 



r(A) = 
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TZ is the socallcd source of convergence: 
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(23) 



(24) 



where R^^, is the Ricci tensor of the metric. is the 
source of shear: 



(25) 



where e" :— Ef + and Cap-yS is the Weyl curvature 
tensor of the metric. As expected, for isotropic metrics, 
like LTB with the observer at the origin, the source of 
shear is vanishing: 
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FIG. 1. Source of convergence as a function of red-shift 
.The sohd curve stands for the FRW model and the dashed 
curve for flat LTB. 



^LTB = 0. 



(26) 



In the case of flat LTB spacetime with corresponding 
solutions (8) the relevant source of convergence is, how- 
ever, non-vanishing: 
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The non-vanishing convergence means that the light 
coming from a source at (r, t) is converged by (27) when 
observed at the origin O. To derive the relevant red-shift 
(z) for each event (t and r) we used the numerical code 
of ( [7]). We have then plotted the convergence as a func- 
tion of increasing z (Fig.l). As it is shown in the picture, 
convergence is increasing with z. This convergence is dif- 
ferent from that of the FRW spacetime. To understand 
this difference one needs to study different observablcs 
such as time delay of images of a source, which goes be- 
yond the scope of this paper, and we will go into its detail 
in a future publication. 



IV. SWISS-CHEESE MODEL OF THE 
UNIVERSE: OBSERVER IN THE CHEESE 

In the previous section we studied a special problem 
where the observer was at the center of inhomogencities. 
To be more realistic we want to study a more general 
problem where the observer is placed somewhere between 
inhomogencities in the universe and observes a distant 
source. The light coming from the source is passing the 
inhomogeneous regions in between and reaches the ob- 
server. For simplicity, we consider two different Swiss- 
cheese models (section II B)) in which the observer and 
the source are in the cheese and the light passes through 
one of the holes: Onion model of Biswas et al. [8] who 
derive a perturbative LTB solution of Einstein equations 



and study the evolution of the density contrast within the 
holes, and the model of Marra et al. [4] who construct a 
non-pcrturbative solution in the holes. In this model the 
universe is completely filled with these holes which form 
a sort of lattice, taking care of the matching conditions 
between LTB and FRW metrics on the boundary of the 
hole. The hole is almost empty except at the bound- 
ary where the matter is concentrated and has an average 
density matching that of FRW density. 
We will trace light rays within these two Swiss cheese 
models coming from a source in the Homogenous Fried- 
man background, the cheese, passing a lens, a hole, and 
finally detected by an observer in the cheese. Deriving 
the source of convergence and shear in this case, is not 
a simple job due to the difficulty of tracing light rays 
in the Swiss cheese which is not as straightforward as 
in the case of the LTB model with the observer at the 
origin. For arbitrary observer we expect a non-vanishing 
shear. The effect of the shear is best studied through 
observables such as radial and tangential arcs [13] of the 
distant galaxies. We will first elaborate on some of the 
basic definitions and then go on to calculate the arcs in 
models just described. 



A. Some basic definitions in lensing 



The surface mass density of the lens is defined as 



m) 



p{r, t)dz, 
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where r = \J z"^ + ^ is the coordinate aligned with the 
line of sight and ^ are the coordinates in the lens plane 
(with the origin on the lens). Convergence of the light 
bundles made by the lens is defined as follows: 
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[13]. 
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Now consider that we are observing a far galaxy and 
the light coming from that galaxy to us is bent due to 
the lensing of the inhomogeneity which exists in their 
path. There are some places in the lens plane where the 
magnification of the lens goes to infinity. Then what we 
see in the lens plane are some radial or tangential arcs. 
The magnification of a lens is defined as follows [13]: 



11 (X 



det{A) ' 



(30) 



where A is the Jacobian of the transformation of the 
source plane coordinates to the lens plane ones. When 
the determinants of A goes to zero we will have the arcs 
in the lens plane (called critical curves). Therefore, to 
derive possible arcs, one should calculate the eigenvalues 



of A first and sec if it can be zero somewhere on the lens 
plane or not. 

The mass inside radius x is described by the dimen- 
sionless function: 



m{x) := 2 / K{^)d^. 
Jo 



(31) 



It can be easily shown that the eigenvalues of the A ma- 
trix be derived using this dimensionless function [13,6]: 
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Ar and At are called radial and tangential critical 
curves, respectively. When A^ goes to zero the arcs of 
an extended source in the lens plane are radial and when 
A( goes to zero the arcs of an extended source in the lens 
plane are tangential. 



B. Radial and tangential arcs I: Onion Swiss cheese 
model 



To study the arcs, we trace light rays coming from a 
far source (in the Homogenous Friedman background - 
the cheese) passing a lens with perturbative LTB metric, 
as described in section(II A), and finally detected by an 
observer in the cheese. Once we have the density profile of 
such a lens, we can immediately derive the place of arcs 
and see if there exists both radial and tangential arcs 
or not [6]. Different density profiles have been studied 
before. For example the NFW profile [14] can produce 
radial and tangential arcs [6] . 

Density profile of a LTB lens with perturbative metric 
(section (II A)) is given by the following relation, derived 
from the Einstein equations in the small u approxima- 
tion: 
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TiG(Mtf 



(33) 



where A(r) := — ^ . 

It seems that the density is decreasing with time, while 
we expect the structure grows as a function of time. But 
the important point here is that the relevant quantity to 
study the structure growth is the density contrast. 

The background spacctime is a FRW spacetime. 
Hence, one can define the density contrast as the devia- 
tion of the holes density {p{r,t)) from the cheese density 
{< p > (t), averaged over the holes) [5]: 



p(r, t)- <p>{t) 
<P>{t) 



(34) 



defining e(r, t) := R2^'^{Mt)'^/'^ A{r), one can derive the 
density contrast: 
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l + e(r,t)- 



(35) 



As Biswas et al. show [5], A{r) has to be bounded 
to get underdensity (voids) and overdensity (structure) 
regions: 

Ajnin ^(^) Ajyiax ^ ^rain (5 <C ^rnax (36) 

As €{r,t) is an increasing function of time, its sign 
is crucial to get the underdensity and overdensity re- 
gions. A negative e(r, t) gives a positive density contrast 
which grows with time (growing structures) and a posi- 
tive e(r, t) gives a negative density contrast growing with 
time (growing voids). Therefore, Amin should be nega- 
tive, corresponding to an overdense region. 

As it is shown in [5], the following E{r) can satisfy the 
above condition (36): 
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where L is a typical length of the large scale structures 
and A\ is the amplitude of the density oscillations. This 
may be seen by looking at R(r,t): 



where A{t) := i?27^Ai (^Mt 
also be written as 



A{t) L . 2 (T^'T 



, (38) 
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The density profile can 



67r(Mt)2 [l + ^(i)sin(2ji:)] 



(39) 



showing the significance of A(t) as the amplitude of os- 
cillations. 

The special characteristic of this Swiss cheese model is 
that each hole has an onion like density profile which at 
the large r's goes to a homogeneous background density. 

Now, consider one of these LTB regions as a lens. The 
plane perpendicular to the line of sight of the source 
(which can be a far galaxy) is called the lens plane. The 
light passes through the FRW region and, close to the 
lens, it is bent due to the density profile of the LTB 
hole using the small u approximation, just as in the 
Schwarzschild case. 

In realistic cases one may forget about the time evolu- 
tion of the lens during the passage of light. This means 
the characteristic time of lens evolution is much greater 
than the light passage time. In the following we will ap- 
ply this approximation. 

We have now all the quantities to calculate A^ and \t 
for the density profile (39). Figures 3-6 show the behav- 
ior of m[x\ Ar, and At as a function of the distance to 
the center of the lens in the lens plane, x. 
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FIG. 2. Surface mass density of the lens as a function of 
the coordinates on the lens plane. 




FIG. 5. Behavior of Xt{x) with respect to x in the Swis- 
s-cheese model. 




FIG. 3. Mass of the lens inside radius x in the sinusoidal 
Swiss-cheese model. 



As expected, m{x) increases as a function of x. Ob- 
viously there is no solution to the equations = and 
Af =0, as can be seen from the Fig 4 and Fig 5. A,. = 
may have a solution in the large 'x' which is by far out 
of the range of our approximation. 

The result of no arcs, neither radial nor tangential, 
means that the onion Swiss cheese model with sinusoidal 
solution is ruled out by observations. 
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FIG. 4. Behavior of \r{x) with respect to x in the Swis- 
s-cheese model. 



C. Radial and tangential arcs II: Swiss cheese model 
of Marra et al 

We assume again that the time of light passage through 
the hole is negligible relative to time evolution of the den- 
sity of the hole in the model of MKMR [4]. To derive the 
solution of Einstein equations within the holes obeying 
the junction conditions, they choose the initial density 
function to have a gaussian profile: 

p{r, U) = Aexp[-(r - rM)^/2a^] + e (r < r^) 
p{r, U) = pFRwiU) {r > rn), (40) 

where e = 0.0025, ru = 0.42, a = r,,/10, vm = .037, A = 
50.59, and pFRwiU) = 25. The hole ends at = .042 
which is equivalent to 350 MPc. This is not a big hole 
but is almost an empty region: the matter density in the 
hole is lO""^ times smaller than in the cheese. 

Applying this initial condition to the Einstein equa- 
tions (in the curved LTB case: E{r) ^ 0), one gets v(r^ t) 
(the peculiar velocity), R{r,t), and p{r,t). Hence, we 
have the density profile of the lens ai t = which, in 
their notation, is the present time (Fig 6). 

We have done the calculation along the same line as 
in the last section, assuming the source and the observer 
far from the hole and in the FRW background (cheese). 
As can be seen from the Figs 7 and 8, this model allows 
both radial and tangential arcs. However, this happens 
at r > rfi which means that the arcs will be observed out 
of the inhomogeneous region. 



V. CONCLUSION AND DISCUSSION 

Independent of how successful the inhomogeneous 
models are in explaining the dark energy problem, grav- 
itational lensing may serve as a criterion to distinguish 
inhomogeneous cosmological models. Different concepts 
developed in the cosmological gravitational lensing tech- 
niques such as shear, convergence, tangential and radial 
arcs, and time delays maybe used to see how tenable 
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FIG. 6. Behavior of 7?(r,t)(y(r,i)) in MKMR model, pecu- 
liar velocity v{r,t) and the density profile p{r,t) with respect 
to r for the curved case at ti — —.8 and t = (present time). 
The straight lines in the R{r, t) diagram are FRW solutions 
[4]. 



these models are. The widely discussed LTB cosmologi- 
cal models, having a vanishing shear as the FRW models, 
differ from FRW ones in the value of convergence which 
may lead to observable effects such as different time de- 
lays of the cosmological images and the large scale lensing 
effects in the CMB. Assuming the observer outside the 
center of symmetry of LTB, one expect a universal shear 
not seen in the FRW models. 

The Swiss cheese models provide us with a density pro- 
file for a " hole" , to be compared with the NFW profile. 
Therefore, the question of tangential and radial arcs may 
lead us to a test of such models, or to a better fixing of 
the model parameters. The onion model predicts neither 
a tangential nor a radial arc. We may therefore rule it out 
even as a toy model to explain dark energy. The MKMR 
Swiss-cheese model [4] do produce both radial and tan- 
gential arcs. The arcs are located in the cheese outside 
the hole near the massive shell. The size of the hole is 
about 350 Mpc, much bigger than familiar structures in 
the universe. Therefore, it is not possible to compare 
this result with real data. It may be possible, however, 
to fix the parameter of the model such that more realistic 
arcs results. If these parameters are compatible with the 
explanation of the dimming of the supernovas is another 
question. We therefore conclude that it is desirable to do 
more research on different aspects of gravitational lens- 
ing effects in inhomogeneous models of the universe. 




FIG. 7. Behavior of Xr{x) with respect to x in Marra et al 
model. 




FIG. 8. Behavior of \t{x) with respect to x in Marra et al 
model. 
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